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Abstract 

In this paper we examine the behavior in temperature of the free energy on 
quantum systems in an arbitrary number of dimensions. We define from the 
free energy a function C of the coupling constants and the temperature, which 
in the regimes where quantum fluctuations dominate, is a monotonically in- 
creasing function of the temperature. We show that at very low temperatures 
the system is controlled by the zero-temperature infrared stable fixed point 
while at intermediate temperatures the behavior is that of the unstable fixed 
point. The C function displays this crossover explicitly. This behavior is rem- 
iniscent of Zamolodchikov's C-theorem of field theories in 1+1 dimensions. 
Our results are obtained through a thermodynamic renormalization group 
approach. We find restrictions on the behavior of the entropy of the system 
for a C-theorem-type behavior to hold. We illustrate our ideas in the con- 
text of a free massive scalar field theory, the one-dimensional quantum Ising 
Model and the quantum Non-linear Sigma Model in two space dimensions. 
In regimes in which the classical fluctuations are important the monotonic 
behavior is absent. 
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1 Introduction 



Our understanding of the behavior of systems near critical points in 1+1 dimen- 
sions has been significantly improved in recent years due to the application of the 
methods of conformal field theory One of the aspects of this method which 
makes it particularly powerful and elegant is that it relates in a simple way the 
microscopic behavior of such systems with thermodynamics. At the critical point 
the most important parameter of the theory is the central charge, c, which is re- 
lated with the quantum fluctuations at zero temperature. The central charge can 
be calculated in many ways. At zero temperature it can be calculated from the 
correlation function of the energy momentum tensor or from the commutation re- 
lations between the generators of the Virasoro algebra 0]. It can be shown that c 
can also be obtained at finite temperatures from the specific heat of the system and, 
therefore, it can be measured in an experiment. 

Conformal Field Theory (CFT) describes the behavior of systems at fixed points. 
Away from fixed points, CFT can also describe the "early" an "late" stages of the 
crossover between fixed points. The behavior at intermediate scales is, in general, 
non-universal and thus depends on the microscopic description of the system. The 
universal behavior of this crossovers is described by the famous Zamolodchikov's 
C-theorem B . This theorem states that it is possible to define a function C of the 
coupling constants of the theory which decreases monotonically under the renor- 
malization group flow, that is, as we change the characteristic scale the function 
C goes monotonically from a fixed point of larger central charge to one at smaller 
central charge. The assumptions present in his proof are the euclidian invariance, 
the existence of a symmetric, conserved, energy momentum tensor and the unitarity 
of the field theory. The physical interpretation of this theorem is related with loss 
of information due to the coarsed graining procedure present in the renormalization 
group equations. Another way to interpret this result asserts that C measures the 
number of fluctuating degrees of freedom. These must decrease when we integrate 
over small scales in the renormalization group flow. 

All the results of the conformal field theory seem to be very difficult to be applied 
for higher dimensions due to the proliferation of components in the energy momen- 
tum tensor. Nevertheless some importants attemps have been made with the use 
of the spectral representation for the correlation function of the energy momentum 
tensor Q, ||. However, the relation between these results and thermodynamics 
remains unclear up to now. 

The approach of this paper is based on the behavior of the thermodynamic func- 
tions under the renormalization group flow. Following this route we can determine 
the macroscopic conditions which are necessary for any system, in any number of 
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dimensions, to have a behavior which is analogue to the one described by the C- 
theorem proposed by Zamolodchikov. We can show that the the monotonic behavior 
implied by the C-theorem can not be derived from the laws of thermodynamics alone 
and, therefore, this behavior is not universally true. As we will show, the sort of be- 
havior needed for a C-theorem to hold is only possible in regimes in which quantum 
fluctuations dominate. In 1+1 dimensions this is always the case because of the triv- 
iallity of one-dimensional classical statistical mechanics. But, in higher dimensions, 
there is a complex quantum-to-classical crossover which spoils the general validity 
of a C-theorem. To put it differently, this is a consequence of the fact that 1+1 
dimensions all systems are at or below their lower critical dimension. This fact im- 
poses stringent constraints on the allowed behavior of the entropy of these systems 
at low temperatures. However, in higher dimensions a variety of complex behaviors 
are allowed and the constraints on the entropy are much weaker. We believe that 
this observation is a hint that there is no general C-theorem in dimensions higher 
than 1+1. 

We have investigated these issues in the context of several systems. We first 
consider the case of a free massive real scalar field in arbitrary dimensions. The re- 
sults also apply to more general bosonic systems as well as to the case of relativistic 
fermions. Here, the crossover parameter is the mass. We find monotonic behavior in 
any number of dimensions. We also discuss two more non-trivial models. The first is 
the quantum Ising model in 1+1 dimensions. We show that the the monotonicity of 
the C-function (calculated without approximations) defined from the themodynam- 
ics is obeyed in the low temperature limit and that near the scaling limit (where 
the model can be described as a field theory of free massive fermions) we found 
all the results expected from the point of view of conformal field theory. At tem- 
peratures comparable or higher than the fermion "bandwidth" , there is a crossover 
to classical behavior and monotinicity is lost (this is the "Debye temperature"). 
Next we discuss the example of the Non-linear Sigma model (the continuum ver- 
sion of the quantum Heisenberg antiferromagnet) in 2+1 dimensions. We show that 
the C-function (calculated from the 1/N expansion) is not monotonic for arbitrary 
temperatures and coupling constant due to a crossover from classical to quantum 
behavior. This crossover can be observed even near the phase transition. Once 
again, in the quantum regime we find a term in the free energy which displays the 
same monotonic behavior expected from a C-theorem. 

In order to understand these results, we observe that temperature plays two dif- 
ferent roles in these problems. At low temperatures, these systems are essentially 
quantum in d + 1 dimensions. Here the temperature mainly affects the geometry in 
which the system lives since it measures the "perimeter" of Euclidean space-time in 
the (imaginary) time direction. As the temperature is raised, this perimeter grows 
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smaller and the temperature drives the system away from criticallity, just like a 
relevant operator. On the other hand, at high temperatures the dimension in the 
imaginary direction effectively shrinks to zero (compared with all the length scales 
in the system) and the system behaves like a classical system in d dimensions. In 
this case the temperature only plays the role of a coupling constant in classical sta- 
tistical mechanics. Thus, at these high temperatures, the C-function of the system 
is not necessarily monotonic. For instance, one does not expect monotonicity if the 
classical system has a phase transition at a non zero temperature. Even if this is 
not the case, such as in classical systems at or below their lower critical dimension, 
non-monotonicity may still occur if the classical system has gapless configurations. 
We will see below that this is the case of the Non-linear Sigma Model in the "classi- 
cal regime". In the case of systems with discrete symetries (e.g., the Ising Model), 
non-monotonicity still occurs near what we will refer to as the classical fixed point, 
but not near the critical zero-temperature fixed point. Thus, we do not expect to 
have a monotonic behavior of any piece of the free energy for all temperatures but 
only in regimes in which the quantum fluctuations are dominant. The strategy of 
this paper is thus to look for the regimes of parameter space ( i. e. coupling constants 
and temperature) in which the physics is dominated by the quantum fluctuations. 

In the next section we discuss the relation between the Renormalization Group 
of quantum systems at finite temperature and Thermodynamics. In section 3 we 
discuss our results in the context of free scalar and fermi fields. In section 4 we apply 
our ideas for the case of the quantum Ising model in 1+1 dimensions. In section 5 
we show how our results work for the case of the Non-linear Sigma model in 2+1 
dimensions. The last section contains our conclusions. 



2 The Renormalization Group and Thermodynam- 
ics 

Suppose we are interested in the study of a system in statistical mechanics which can 
be described in a lattice with some characteristic length a (UV cut-off) and which 
has a set of dimensionful coupling constants g. By naive dimensional analysis we 
know that the bulk free energy of this system in thermal equilibrium at temperature 
T (P = 1/T, H — ks — 1) in (d + l)-dimensions, in the thermodynamic limit, can 
be written as 

mg , a) ^ EoM .m^iA (2 ,) 

where E is the zero temperature energy, V is the volume (V — > oo, but N/V is 
finite, where N is the number of particles), v is a characteristic velocity in the system 
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(for instance, the velocity of the quasi particles), n(d) is a positive real number which 
depends only on the dimensionality of the system (see next section for details) and 
C is a dimensionless function of its arguments. 

We want to comment here that in (2.1) we are implicitly assuming that the 
ultra-violet divergences in the field theory can be subtracted in the first term on 
the r.h.s. of (2.1) since these divergences only occur at zero temperature, that is, 
since C is a finite size effect (in terms of temperature) it is insensitive to divergences 
which occur in small scales |J. 

Here we are interested in the limit when a — > 0, that is, in the limit where the 
system can be described by a field theory (we do not intend to discuss here how this 
can be achieved, we just assume that is possible). In this limit the "bare" coupling 
constants of the original theory are replaced by renormalized dimensionless coupling 
constants a and since C is dimensionless we should write 



where A has dimensions of inverse of length. 

Since C is an observable it does not change as the scale is changed and therefore 
it obeys a renormalization group (RG) equation (t = (/?A) _1 ) 



is the beta function of the system. Observe that the beta function can also be 
obtained from (2.3). 

The solution of (2.3) is obtained with the introduction of the correlation length 
of the system, £(ck, t), defined by its renormalization group equation |7| 




(2.2) 




(2.3) 



where 




(2.4) 




(2.5) 



It is trivial to show that 




(2.6) 



Define now a dimensionless correlation length by 



7 = A£ 



(2.7) 



5 



and therefore from (2.6) 

C = C(t-f(a,t)) (2.8) 

and from (2.5) we have 

^ -*«>£■ <»> 

It is also important to define the running coupling constant, a(t), as the solution 
of the equation (2.4) 

t^= -P(<*(t)) (2.10) 

with the initial condition 

a(t = t ) = a (2.11) 

where to is an arbitrary renormalization point. 

It is well known that at a critical point (where the conformal field theory can be 
applied 0) the function C as defined in (2.1) is a constant and, in one dimension, 
is related with the Virasoro algebra (it is called the central charge of the theory). 

Zamolodchikov || proposed some years ago a function Cz(R,g) (where R is a 
length scale in the 1+1 dimensional space) which is monotonically decreasing under 
the RG flow, that is, 

< 2 - 12 > 

and which, at the critical point, is the central charge. His proof is based on the 
euclidian invariance, the existence of a symmetric, conserved, energy-momentum 
tensor and unitarity of the field theory. 

Suppose we can actually calculate the spectrum for the system under considera- 
tion. Of course the free energy can always be written in the form (2.1) but, in this 
case, since C is dimensionless, 

C = C{<f>(g)0) (2-13) 

where <fi(g) is a characteristic energy scale which depends on the bare coupling 
constants of the theory and possibly on the temperature. In general we are talking 
about a Hamiltonian with the form 

H = H + gH 1 (2.14) 

where Hq is the fixed point Hamiltonian and Hi is a "pertubation" which drives the 
system away from criticallity. In this case when g vanishes C goes to the central 
charge at the fixed point, 

Lim g ^oC(<j)(g)f3) = Co = constant. (2-15) 
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Since in thermodynamics, for each system, the coupling constants are fixed, the 
only energy scale free to be changed is the temperature. Therefore we proposed an 
analogue of the Zamolodchikov's equation (2.12) which is valid in any number of 
dimensions 

Kf)/° (2 - i6) 

Observe that we are using the fact that in finite temperature field theory we 
have to impose periodic (or anti-periodic) boundary conditions in the imaginary 
time direction. These boundary conditions permit us to interpret a quantum system 
in d + 1 dimensions as a classical system in D(= d+ 1) dimensions. The manifold 
under consideration has the form of a cilinder since the other d dimensions are free. 
The length of the cilinder in the imaginary time direction is (3 and, therefore, (2.12) 
and (2.16) represent the same change in the scale, although, we could say, (2.12) is 
a isotropic change while (2.16) is anisotropic. Observe that at high temperatures 
the length of the cilinder goes to zero and therefore the field theory is restricted to 
d dimensions. 

From the point of view of the renormalization group approach we see that if we 
follow the flow using the running coupling constant defined by equation (2.10) the 
function C is stationary as required. However if we fix the renormalized coupling 
constant (as required in Zamolodchikov's theorem) and vary the dimensionless tem- 
perature, t, we would get (2.16). Using the result (2.8) and the RG equation (2.9) 
we find, 

dC^J^d^ dC d 1 ~ 

dt d(t 7 ) dt d(t-f)da P{) { ] 

This expression will be useful later. 

Consider a system initially at temperature zero in a volume V, with ground state 
energy E and entropy So. Suppose we perform a thermodynamic transformation at 
a constant volume which increases the temperature of the system up to a tempera- 
ture T. At this point the system, in thermal equilibrium, has an internal energy U 
and an entropy S. We will show that if, 

XTT 

(2-18) 
where 

5U=U-E (2.19) 

and 

SS = S - S (2.20) 
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then (2.16) follows. We will show at the end of this section that (2.18) can not be 
obtained from the laws of thermodynamics. Therefore the irreversibility of the RG 
flow expressed in (2.16) is a consequence of (2.18). 

Notice that (2.18) has a very nice interpretation in terms of the number of states. 
Since in the thermodynamical limit the entropy of a system can be writen as || 

S(U) = ]n(fl(tf)) (2.21) 

where Q(E) is the thermodynamic number of states defined as the number of states 
with energy between E and E + dE, we can reexpress (2.18) as 

{l(E )e^^ > n(U) (2.22) 

which means that the number of states with energy between U(T) and U(T) +dU(T) 
is bounded from above. This is a highly non trivial statement which, as we shall see 
below, is not required by the laws of thermodynamics alone. 

Before we start to prove that (2.16) is a consequence of (2.18) let us remember 
that the continuum limit a — > is very important in Zamolodchikov's proof. There, 
he assumed that, 

R > a (2.23) 

In this paper we also assume the counterpart of (2.23) in terms of temperature, 
namely, 

T < 6 (2.24) 

where 9 is a characteristic cut-off temperature (for instance, the Debye temperature, 
Od). 

Using equation (2.1) the thermodynamical quantities, pressure (P), entropy (S) 
and thermal energy (U) can be easily calculated, 

n{d)C 



n{d)V A _ Q dC s 
<d)V ( dC* 



S = ^w[(d+l)C-P^) (2.26) 



U = Eo + ^[dC-(3^) (2.27) 



Notice that as the temperature goes to zero the free energy and the internal 
energy goes to the ground state energy. This fact means that C and dC/d(3 go to 
zero in this limit. From the third law of thermodymanics we know that 

Lim T ^ S(T) = S (2.28) 
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So in order to our definition of C be consistent we have to set 



S = (2.29) 

which is allowed from the thermodynamic point of view. 
From (2.26) and (2.27) we can write 

r)H i) d Pi d 

^ = ^k (dS - [i+l)m) (2 ' 30) 

therefore if (2.18) is obeyed (2.16) follows and we get the same monotonicity prop- 
ertie as in the original C-theorem. 

We can also rewrite (2.30) in an useful form. Recall that the definition of the 
specific heat at constant volume is 

C V (T) = T^ = ^- (2.31) 
which, when integrated using (2.29), gives 

S(T) = ^jp-dT' (2.32) 

and 

5U(T) = f C v {T')dT'. (2.33) 
Jo 

Substituting (2.32) and (2.33) in (2.30) we find 

Expression (2.34) is very suitable for experimental test and can be used to classify 
the systems which obey assumption (2.18). Morover, at the fixed point we expect 
that dC/dp = and therefore the specific heat must behave like T d which is known 
to be true in one dimension. 

Now we can explain why condition (2.24) is important for our results. The 
pertubation term in (2.14) breaks the scale invariance of the fixed point hamiltonian 
and it drives the system into the basin of attraction of a more stable fixed point. 
Once temperature is turned on, this operator generates flows away from the fixed 
point but which are "cutoff" once the correlation length £ is of the order of the 
inverse temperature. Thus, at temperatures higher than the system behaves 
as if it were still at the infrared unstable fixed point but, at lower temperatures 
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it crosses over to the infrared stable one. It is this crossover which is described 
by Zamolodchikov's theorem ||. Morover, at high temperatures the specific heat 
must go to its classical value which is a constant given by the equipartition theorem 
(Dulong- Petit law), therefore, the entropy diverges logarithmically and we do not 
expect that assumption (2.18) is obeyed. However, if the spectrum has a width which 
is much greater than the energy scale 4>(g) we could imagine a situation where the 
temperature is much greater than <j)(g) but much smaller than the bandwidth. In 
this case, since C is only a function of (fi(g)/3, the limit as the temperature goes to 
infinity is identical (in terms of RG) to the limit as g — > and therefore, by (2.15), we 
see that the C is the central charge in this point. This is a stable high temperature 
fixed point. Nevertheless, for temperatures larger than the bandwidth we expect 
that the behaviour of the system is classical and (2.18) will not be obeyed. We will 
see later that this behaviour can be seen in the Ising model in 1+1 dimensions. 

For systems with a phase transition (when also depends on the temperature) we 
know that near the critical temperature the correlation length of the system increases 
exponentially and therefore we do not expect that the thermodynamic number of 
states be bounded from above. This behaviour near the critical temperature is 
essentially classical and therefore the monotonicity of C should be violated. We will 
show that this is what happen in the Non-linear Sigma Model. 

We want to point out that in the assumption that we can describe the system by 
a quantum field theory we are assumig that there is a zero temperature fixed point 
which is unstable under the increase of the temperature. This is clearly expected in 
any quantum system since we always have the crossover from the quantum to the 
classical system as the temperature increases. Therefore, (2.16) shows that the RG's 
flow goes from the zero temperature fixed point to a fixed point at finite temperature 
( in particular to a high temperature fixed point). 

In order to finish this section we will show that we can not obtain (2.18) from 
the laws of thermodynamics. 

Suppose we have initially the system at zero temperature and we put it in contact 
with a thermal reservoir at temperature T and wait for the thermal equilibrium. 
Suppose that this system is isolated from the rest of the universe. The variation of 
the entropy of the reservoir is —Q/T where Q is the heat given by the reservoir to 
the system. Using the second law we get 




(2.35) 



Using the first law we have 



5U = Q 



(2.36) 
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therefore 

S > jr (2-37) 

We see that the laws of thermodynamics can not lead to (2.18), that is, they show 
that the increase of the entropy is bounded from below. Furthemore, notice that 
from (2.37), (2.26) and (2.27) we obtain that C is always positive as required. 



3 Free Massive Field Theories 

We will now discuss the physical meaning of the results of the last section in the 
context of specific systems. We will first consider the case of a massive scalar field 
at non-zero temperature and exhibit the crossover. We will do that by showing that 
the restriction (2.18) is always obeyed. 

The approach is very straightfoward. The free energy for a system with one- 
particle excitation spectrum E K where k is a set of quantum numbers is given by 



where r is +1 (-1) for bosons (fermions). 
Comparing (2.1) and (3.1) we get 



F = -Bo + ^ 5> (l - re--" 5 -) (3.i; 



C = -g£rS>(l-™-~) (3.2) 



Suppose we have the relativistic dispersion relation 

E k = Vm 2 v 4 + v 2 k 2 (3.3) 

where k is the wavenumber (in units of h) which is defined by periodic boundary 
conditions. In the thermodynamic limit we can replace the sum in (3.2) by an 
integral as usual. We will also assume v — 1 and rewrite (3.2) as 



C = --J——T I d d k In f 1 - re-^ k A (3.4) 

n ^ 9-ttV \ J V ' 



— * 

since the spectrum only depends on the absolute value of k, we have 
where T(x) is the Gamma function. 
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Using (3.3), changing variables in the integral (3.5), expanding the logarithm in 
powers and integrating, we easily get 



4±i oo l+l 



c W) = A L m E tet* *ti c ( 3 - 6 ) 

n(d)2 2 it 2 i = i I 2 2 

where K u (x) is the Bessel function of imaginary argument of order v. 

Observe that, as expected, C is only function of one variable, namely m(3, and, as 
explained before, we find that the correlation length is given by m~ l . Furthermore, 
since a = A _1 m ( 7 = a" 1 ), the beta function is trivial, namely, /3(a) = —a, and 
the running coupling constant as defined in (2.10) is just a(t) = at/to- 

We can also obtain the special cases. In the massless limit, when kinetic energy 
of the system (given by m) is much smaller than thermal energy (given by (3~ l ) 
that is, m(3 « 1 (£ >> 1/T), we can approximate the Bessel function for small 



argument [10 and write 



n(d)7i 2 l=1 1 

which is finite. This is the high temperature fixed point. 

From (3.6) and using well know relations between Bessel functions [10] we find 
that the rate of change of C with the temperature is given by, 



dC 
~df3 




T K^ 1 (lm(3) (3.i 



which is always negative in any number of dimensions. We could also arrive to this 
result using equation (2.17). 

However the normalization factor n(d) is still arbitrary. For relativistic systems 
in 1+1 dimensions, this factor, at a fixed point, equals to the central charge of the 
Virasoro algebra which can be thought as a measure of the number of fluctuat- 



ing fields. Its well known ]TT| that in 1+1 dimensions relativistic massless bosons 
have C(0) = 1 and relativistic massless spinless fermions (Majorana fermions) have 
C(0) = 1/2. Hence, 

n(l) = |. (3.9) 

In higher dimensions we do not know if a generalization of the notion of central 
charge exists. Still, we may choose the normalization factor in such a way that, 
for free fields, it counts the number of degrees of freedom. This choice leads to the 
definition 

T (d±i) 

<d) = -^((d+l) (3.10) 

7T 2 
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for bosons and 

r(*£) (2 d -i) 

n(d) = V ^ + 1) (3.H) 

7T 2 2 

for fermions, where C(a;) is the Riemman's Zeta function. Observe that spinless 
fermions and bosons only have the same n(d) at d — 1. Anyway this arbitrariness in 
n(d) does not affect our discution about the mononicity properties of the C function. 

Let us remark that the behaviour of the function C is non analitic as a function 
of m, that is, if m is strictly zero the function C is a constant given by (3.7) in all 
ranges of temperatures. If m finite, even very small, the behaviour of C is completely 
different (see Fig.l). Indeed, in the massive limit, m/3 » 1 (£ << 1/T), or low 
temperatures, the Bessel function in (3.6) vanishes exponentially and C is zero. 

Finally, for systems like Fermi liquids we know that independent of the dimen- 
sionality the specific heat is linear in the temperature []E|. If we substitute this 
result in (2.34) we will see that dC/d(3 is negative in two or three dimensions but it 
is positive in one dimension. This result is clear from the point of view of relativistic 
particles since in one dimension we can always expand the energy around the Fermi 
surface as in (3.3). 



4 The Ising Model 



Here we will consider the Ising model which can be described by the following 
Hamiltonian 



N N 



H = -vY,sl-JY,s 1 n s 1 n+1 (4.i; 



n=l n=l 

■th 



where S 3 n is the j (j = 1,2,3) projection of the spin operator in the site n (n = 
1,...,N). We will assume periodic boundary conditions, that is, 

Sh+i = Sl (4.2) 

and neglect any surface effect which might occur. 

This model was solved exactly [Tj| using the Jordan- Wigner transformation. 
Since the procedure is somewhat standard we will only quote the results (see [|13[ 
for details). 

It is possible to show that the problem reduces to a problem of fermions with 
dispersion relation given by 



E k = TJl + g 2 - 2gcos(k) (4.3) 
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where 

27rm 

* = — < 4 - 4 > 



with m = 0, ±1, ±2, ±y (iV — > oo). The constant £> is defined as 

I = ± (4.5) 

and it controls the physics of the problem. Observe that when g = 1 the dispersion 
relation is simply 

E k = 2Tsm(^j (4.6) 

and the fermions are massless. The velocity of the excitations can be easily computed 
to be T near k = 0. At this point the system is scale invariant and it represents the 
phase transition. 

The free energy is obtained straightfowardly JTJ 

F = E -— dk In (l + e~P E A (4.7) 

p7T JO V ' 

where 

E = -— dkE k (4.8) 

Z7T JO 

is the ground state energy. 

From (2.1), (3.9) and (4.7) we get 

C=^E r dk In (l + e-^) (4.9) 

TT JO V ' 

and although we do not have an analytic form for (4.9) we can study its limits. 
Observe that at the critical point, g = 1 (J = 2r), we can rewrite (4.9) as 

c = QJY ~ (-If r/ 2 ^ e _ 2 ^rsin( fc / 2 ) (410) 
7T 2 ^ n ^0 

and therefore at the critical point C depends on the temperature. This occurs 
because we know that the continuum limit can only be taken at low temperatures. 
Recall that k is the lattice momentum. For a system with lattice constant a, the 
momentum q in laboratory units is q = k/a. The continuum limit is the limit of 
a — > 0, keeping the physical dimensionless temperature t -1 = Ta/T fixed. In this 
limit, (4.10) becomes 



limC = lim 4- £ — f* dqe- 2 ?-^) = - (4.11) 

a ^0 a^O lT 2 t^ 71 Jo 2 



?1=1 
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as expected from the results of conformal field theory. 

We can also understand this result from the point of view of section 3. Observe 
that at low temperatures and close to the phase transition we can expand (4.3) 
around k = as 



E k = Ty/(l-Q)* + Qk* (4.12) 

which has the same form as (3.3) where the mass term (or the inverse of the corre- 
lation length) is proportional to 1 — g which vanishes at the critical point. If now 
we let k varies between and infinity we recover the result of the former section. 

Observe that here the fact that the energy cuttoff ( i.e. the bandwidth of the 
spectrum) is finite is very important. We can not excite particles in the system with 
an energy larger than the bandwidth and therefore we have a temperature cut-off 
exactly as discussed in second section for the case of constraint (2.24). This fact 
will result that the C function defined in (4.9) will decrease with temperature for 
temperatures larger than the bandwidth. If we take the limit that (3Y « 1 in (4.10) 
we easily see that C goes to zero in this limit. 

We will pay attention to two limits, the limit where the we have unbroken sym- 
metry (the vaccum is invariant under the symmetry of the Hamiltonian) and the 



limit of broken symmetry | 15[| . When g « 1 (J <<T) the vaccum is an eigenstate 



of S 3 where all spins are pointing up. From (4.9) is easy to see that 

C = ^£ In (l + e^ r ) (4.13) 



and get a behaviour where C increases monotonically in the low temperature re- 
gion and when reachs temperatures of order of the bandwidth, namely T, C goes 
monotonically to zero (see Fig. 2). 

When g >> 1 (J >> T) the vaccum is an eigenstate of S 1 and it is degenerate 
(all the spins can be up or down), that is, the symmetry is broken. From (4.9) is 
trivial to show that 

(7 = ^ In (1 + e"¥) (4.14) 

and therefore we see that C increases monotonically from T = (C = 0) to tem- 
peratures of order J and then decreases monotonically to zero again. We see that 
the function C in the low temperature regime increases monotonically with the tem- 
perature as expected by Zamolodchikov's theorem and the cut-off temperature, as 
explained in (2.24), is J. This result can also be obtained from the free energy of 
a classical one-dimensional Ising model. It is well known that this trivial classical 
model has an infrared unstable fixed point at T = 0. We call this the classical zero 
temperature fixed point. In this classical problem it is obvious that the free energy 
is not a monotonic function of the temperature even though the correlation length 
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is indeed monotonic. Monotonicity is only present near the fixed point, i.e. near 
T = 0. 

Anyway, even away from the critical point the monotonicity properties of the 
C function are consistent with Zamolodchikov's theorem if we are in the region 
where the temperature is low compared with the temperature cut-off of the problem, 
exactly as considered in (2.24). Again the behaviour of the C function is non analitic 
in g, this is a clear sign of problems in one dimension. 



5 The Non-Linear Sigma Model 

The non-linear sigma model is the described by the following Lagrangean density 



1 

2£ 



L = — (9^n(S,t)f (5.1) 



with jj, — 1, d+ 1 and the constraint 



(n) 2 = 1 (5.2) 

where n in a vector with N components and g is the coupling constant. 

It can be shown that the Heisenberg antiferromagnet in the continuum limit can 
be described by the non-linear sigma model ( see and references therein). This 



model has been studied in the context of the High- Temperature materials with the 



use the renormalization group method ||17|| . Here we will use the expansion 1/JV in 



order to study its properties and we will see that our results are equivalent to those 
found in |L7| and |L8| . 



The thermodynamical partition function is written as 

Sojn] 



n 2 - 1J e — (5.3) 
where r is the imaginary time and 

ff3 r+OO 

S Q = dr dx d L(n(x,r)) (5.4) 

JO J-oo 

is the euclidean action with — 1/T is the inverse of the temperature. 

Observe that in the low temperature limit we can extend the integral (5.4) over 
the whole imaginary time axis and we find a non-linear sigma model in d+ 1 dimen- 
sions with the coupling constant g. In the high temperature limit the upper limit 
in the integral in (5.4) is small and we can rewrite the action as the action for the 
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classical non-linear sigma model in d dimensions with a coupling constant given by 
We now decompose the field n in its components as 



n — (a, 7r 



(5.5) 



where tt is a vector with N — 1 components. Substituting (5.5) in (5.3) changing 
variables in the functional integral as if = yfg<j) and tracing over <p we get 



Z = j Da(x,r) J DX(x, r)e 



S[<t,A] 



where we have introduced a representation for the delta function and 



S 



dr I d d x - 

OO 2 



{d^af + A (a 2 - 1 



+ 



N -1 



tr In (-& + X 



(5.6) 



(5.7) 



where tr is the trace. 

Now rescale the parameters as 



§=(N-l)g 



a = J(N-l)g a 



and rewrite (5.6) as 



Z = J D~o J DX e -(^-i)^/P,A] 
where the effective action is given by 



dr I d x - 

-oo 2 



(d,a) 2 + Xa 2 -^ 



+ tr In ( -d 2 + X 



(5.8) 
(5.9) 

(5.10) 
(5.11) 



In the limit as iV — ► oo the solution of (5.10) can be obtained by the saddle point 
equations 

(-dl + X)a = (5.12) 

& 2 = — — G\ (x, t; x, t) (5.13) 
9 

and 

(-81 + A) G x (x, t; x>, r') = S d (x - x') S (r - r') (5.14) 
Assuming that a and A are constants we can replace the equations above by 

Xa = (5.15) 
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and 

^'-k/^I^ (5-16) 



P^J (27r)^ 2 +^ + A 

where we have used the periodic boundary conditions for the fields in the imaginary 
time direction 

${x,t) = $(x,t + P) (5.17) 
and therefore w„ = with n = 0, ±1, ±2 



We can proceed further and perform the sum in (5.16) |10| 

2 Jo VWT\ 

where S^ 1 = 2 d - 1 7rlr(f). 

From (5.15) is easy to see that we have two phases in the problem, one ordered 
phase (A = and a / 0) and one disordered phase (A ^ and a = 0). In the 
ordered phase we have N — 1 Goldstone bosons and one massive boson and in the 
disordered phase the full symmetry of the system is recovered and the modes are 
massive. 

Notice that expression (5.18) contains a ultra-violet divergence (k — > oo) at zero 
temperature in one and two dimensions which means that we have to renormalize 
the theory. In order to do so define the renormalization group transformations 

o = ^M (5.19) 

A = |l A 2 m 2 (5.20) 

g = Z^ l - d a (5.21) 

T = At (5.22) 

U = gT = ZxA 2 - d u (5.23) 

where 

it = cat (5.24) 

Here M, m, a, t and u are dimensionless quantities and A is a scale with dimensions 
of inverse of length (A = 1/a, where a is the lattice spacing). Observe that the 
temperature T is not renormalized since it fixes the length of the manifold (the 
cilinder). In order to achieve the renormalization in two dimensions we have to 
introduce a new variable U, which is the high temperature coupling constant as 
explained before. 
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It is easy to see that the theory is renormalized in any number of dimensions if 
we choose 

Z Q = Z t (5.25) 

and 

i , aK d r oath (£V£*+r 



1 + ^/ dxx" "\ ' (5.26) 

Z 2 Jo \Jx 2 + 1 

From (5.26) is possible to calculate the /3-functions of the model. Using (5.21), 
(5.22) and (5.23) we have 



#* = A '!x) = (d-l)a-an(a,u) (5.27) 



{ du\ 

/3„ = A — = (d-2)u-u Sl(a,u) (5.2? 



dA 



g,T 



and 



where 



A = A(*) =-t (5.29) 



(dlnZ Q \ ~ (dlnZ \ ~ fdlnZ \ 

We can now solve the system formed by (5.27) and (5.28) using expression (5.26). 
Although we can do it in general we will choose d = 2 since all integrals can be done 
explicitly. For d = 2 we find (the zero temperature integrals are done as in [pij 
the context of e expansion) 

a 2 , / a \ 



in 



A = Q __ coth( _j (5.31) 
and 

* ^ -£ coth (£) (5 ' 32) 

which agrees with the one loop result found in |DJ. Furthermore, we can solve (5.18) 
explicitly and one finds (use definition (5.24)) 

/ 27r(A/ 2 a-l) / 1 \ \ 

m = 2t arcsinh I e at sinh I — J J (5.33) 

which, together with 

Mm = (5.34) 
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completes the solution of the saddle point equations. 

Notice that the fixed points of the model (where the the f3- functions vanish) are 
u = 0, a = and u = 0, a = 4ir (or simply, (0,0) and (0,47r)). (0,0) is a stable 
fixed point and (0, 4n) is unstable. 

Observe that at zero temperature in the ordered phase (m = 0, M / 0) we get 

M 2 = 1 - — (5.35) 

where 

a c = 4vr (5.36) 

is the critical dimensionless coupling constant which was expected from the above 
analisys of the /3-functions. Observe that the expression (5.35) is only valid as 
a < a c . In the desordered phase (m ^ 0, M = 0, a > a c ) we get 

m = 1 - — . (5.37) 

a 

We can also solve exactly the equation (2.9) for the running coupling constant 



a(t) = at \t + ^- ln( ^T) 1 1 " (5 ' 38) 



We can easily see that the dimensionless correlation length as defined by (2.7) is 
mT x (use (5.33) in (2.9)). Another possible way to achieve this result is to expand 
the action (5.11) around the saddle point solutions and recall the definition of the 
free energy 

F = ~\nZ (5.39) 

we get that the finite term (which does not need to be regularized) in the leading 
order term in 1 /N is given by 

^7 - - S 4S^ K /2 ) * 1 4r {^) (5-40) 

The free energy only depends on the mass of the N — 1 modes, which is represented 
by A. Notice the resemblance with expression (3.6), the main difference is that here 
the correlation length depends on the temperature. Morover, comparing (5.40) with 
(2.1) we find the C function, 

S d 2^T(l) (m(a,t)\^ °° 1 (m(a,t)\ 



n 



(d) 7rV2 
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Comparing (5.41) with (2.8) we conclude that 7 = m~ l , as expected. 

Now we want to know how C changes with the temperature. From (3.8) is easy 
to see that 

and from (5.33) one finds 

Ta £ °' (5 ' 43 > 
Using our result (2.17) we arrive to the conclusion that the sign of the derivative in 
(2.17) equals to minus the sign of the beta function. 

We can clearly distinguish two regions in the parameter space. One of them has 
a positive beta function and is called |T7| the classical renormalized region because 



the correlation length diverges exponentially near T = as it does in the classical 
version of the theory with the bare coupling constant replaced by the renormalized 
one (see (5.33)); the second region has a negative beta function and is called quantum 
region. The border between these regions is given by the zeros of (5 a , namely, 

a = a c tanh (J-^j (5.44) 

The meaning of this border line is clear, we see that at low temperatures the corre- 
lation length at this line equals the thermal wavelength (which in our units is 1/T). 
Therefore, in the classical region the correlation length is always greater than the 
thermal wavelength which gives the classical character of the region, in the quantum 
region the correlation length is smaller than the thermal wavelength, as expected. 

We finally conclude that if we start the flow from the classical region our C 
function will decrease with increase the temperature (since /3 a is positive) as far 
as it finds the border line (5.44) and then, in the quantum region, it will start to 
increase monotonically with the increase of the temperature (see Fig.3). 

We see that the non monotonicity of the C is again related with the crossover 
between quantum and classical behaviour. In one dimension the crossover is simply 
dictated by the temperature in a trivial way, that is, any lattice leads to a bandwidth 
in the spectrum of excitations which is a natural temperature cut-off. This behaviour 
is implicit in the original assumptions of the Zamolodchikov's theorem. In higher 
dimensions the existence of phase transitions (even in the presence of a infinite 
bandwidth as in the case of the Non-linear Sigma model) produces a region in 
the parameter space where the behaviour is clearly classical, that is, where the 
correlation length is much greater than the thermal wavelength. 

Now we can understand why the Zamolodchikov's theorem is so powerful in 
one dimension in the continuum limit. In d = 1 we have only a fixed point at 
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(0, 0) and therefore there is no classical region in the parameter space. This is also 
true for systems with relativistic dispersion relation and although we can arbitrarily 
choose a region in the parameter space where the correlation length (independent 
of temperature) is smaller than the thermal wavelength we do not have a phase 
transition, that is, the correlation length does not blow up dinamically. 

6 Conclusions 

We have shown in this article that, in regions of the parameter space (coupling 
constants and temperature) where quantum fluctuations dominate, it is possible to 
define a C function (in any number of dimensions) which has the same properties 
as the C function defined in Zamolodchikov's C theorem. We showed that in the 
regions where the system behaves classically (where the correlation length is greater 
than the thermal wavelength) monotonic properties are not to be expected. 

We illustrate these issues in the context of several models. We showed that free 
massive field theories have the mononicity properties proposed in Zamolodchikov's 
theorem in any number of dimensions. This result is due to the lack of a classical 
region in the parameter space of the theory, exactly as in one-dimensional field 
theories. Applying our ideas for the case of the Ising model in 1+1 dimensions we 
showed that the Zamolodchikov's theorem applies in the region of low temperatures 
(as expected) and the fact that we have a bandwidth (or temperature cut-off) implies 
a crossover between classical and quantum behavior. 

In the case of the Non-linear Sigma model in 2+1 dimensions we found that in 
regions where the correlation length is greater than the thermal wavelength (that is, 
close to a phase transition) the thermodynamic functions exhibit classical behavior 
and consequently the function C is not a monotonic function of the temperature. 
Nevertheless, in regions where the system is essentially quantum mechanical our 
function C does increase monotonically with the temperature, a result similar to 
the Zamolodchikov's C-theorem. In this way we see that C can be used to describe 
the crossover between classical and quantum behaviour. 

Since our approach is thermodynamical, it is possible to extract these behavior 
from experiments which measure the specific heat of the system under consideration. 
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